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1 Introduction. 



Dimension two condensates have received great attention in recent years. These condensates 
might play an important role for the infrared dynamics of Euclidean Yang-Mills theories, as 
supported by the considerable amount of results obtained through theoretical and phenomeno- 
logical studies as well as from lattice simulations [IllllSlll|51ElI71IHliniIiniIIIlIIlII31IIlII3 

HSlinillHllinillllEIlESESlEllESlESl- 

For instance, the gluon condensate (^JJA^) has been largely investigated in the Landau gauge. 
As pointed out in j4j, this condensate enters the operator product expansion (OPE) of the gluon 
propagator. Moreover, a combined OPE and lattice analysis has shown that this condensate 
can account for the l/Q^ corrections which have been reported jl8l 1191 QUI 1211 UM I25j in the 
running of the coupling constant and in the gluon correlation functions. An effective potential 
for (^AI^A^j^ has been obtained and evaluated in analytic form at two loop in {I^ iIT H IT U IT ^ [TK] . 
showing that a nonvanishing value of (^A'^A'Jj^ is favoured as it lowers the vacuum energy. As a 
consequence, a dynamical gluon mass is generated. We also recall that, in the Landau gauge, 
the operator A^j^A^ is BRST invariant on-shell, a property which has allowed for an all orders 
proof of its multiplicative renormalizability. Its anomalous dimension is not an independent 
parameter, being expressed as a combination of the gauge /3— function and of the anomalous 
dimension, 7^1, of the gauge field A'^ jjjj. This relation was conjectured and explicitly verified 
up to three-loop order in j28j . 

The dimension two operator A'^A'^ has been proven to be multiplicatively renormalizable to 
all orders in the more general class of linear covariant gauges j29j . An effective potential for the 
condensate (^^A^) in linear covariant gauges has been evaluated in ^31) providing evidence for 
a nonvanishing value (^JJA^) in these gauges. 

A renormalizable mass dimension two operator can be introduced in other covariant renormaliz- 
able gauges, such as the Curci-Ferrari and the maximal Abelian gauge. In the Curci- Ferrari gauge 
the generalized gluon-ghost operator (^^JJA^ -|- ac^c") is BRST invariant on-shell, displaying 
multiplicative renormalizability to all orders [HOI- The fields c", c" stand for the Faddeev- 
Popov ghosts, while a denotes the gauge parameter. Evidence for a nonvanishing condensate 
(^A^^^ -|- ac'^c'*) have been provided in J2j. Note that in the limit a ^ 0, corresponding 
to the Landau gauge, the operator (^A^^^ -|- ac"'c°') reduces to A mixed gluon-ghost 

operator, namely [^A'^A'^ + ac'^c'^), can be introduced also in the maximal Abelian gauge 
ISlini^^- Here the color index A runs over the A^(A^ — 1) off-diagonal generators of the gauge 
group SU{N), A = 1, A^(A^ — !)• As in the case of the Curci-Ferrari gauge, this operator is 
BRST invariant on-shell, being multiplicatively renormalizable to all orders [51 lUl ITU I3UI I31j . 
Analytic evidence for a nonvanishing condensate (^^A'l^A^ + ac'^c^) in the maximal Abelian 
gauge can be found in We underline that a nonvanishing condensate (^^A'j^A^ + ac'^c^'^ 
gives rise to the dynamical mass generation for off-diagonal gluons, a result of great relevance 
for the so-called Abelian dominance, supporting the dual superconductivity picture for color 
confinement. An off-diagonal gluon mass has also been reported in lattice simulations |321 133j . 

Studies of the influence of these condensates on the gluon and ghost propagators when the non- 
perturbative effects of the Gribov copies are taken into account can be found in j34[ I35 | l36 l I37j . 
The output of these analysis is an infrared suppression of the components of the gluon propagator 
in the aforementioned gauges, a feature in agreement with the results available from lattice and 
Schwinger-Dyson studies [^113311311111121113111211313211331110111311311312011211^ 
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Certainly, many aspects related to the dimension two condensates deserve a better understand- 
ing. This is the case, for example, of the gauge invariance, a central issue in order to give a 
precise physical meaning to these condensates. A recent study of this topic has been given 
in |nSl[M[IS3) where a set of conditions which should ensure the independence of the conden- 
sate (A^^^) from the gauge parameter in the class of linear covariant gauges has been proposed. 

In this work we pursue the study on the aspects of the gauge invariance of the dimension 
two condensates. Our aim here is that of discussing the possibility of introducing a suitable 
colorless dimension two operator 0{A) which preserves gauge invariance 

60{A) = 0, 

= -Dfu;', (1.1) 

where Df!^ is the covariant derivative 

Dab^gabQ^_gfabcj^c^_ (1.2) 

This is a difficult task, due to the lack of a local gauge invariant mass term built up with gauge 
fields only. This problem could be overcome by looking at nonlocal operators. However, even 
if we allow for nonlocal operators, we cannot give up of the requirement that a consistent com- 
putational framework, allowing to carry out higher loop calculations, has to be at our disposal. 
This is a strong requirement which, in practice, deeply constrains the type of nonlocality allowed 
for the dimension two operator. As a suitable proposition in order to obtain such a consistent 
framework, we could demand that the action to which the nonlocal gauge invariant operator 
0{A) is coupled, should have the property of being made local by the introduction of a suitable 
set of additional fields. 



(I.) Therefore, denoting by So the term which accounts for the introduction in the Yang- Mills 
action, Sym, of the operator 0{A) in its localized form, we require that So is gauge 
invariant. 

(II.) Also, on physical grounds, we demand that the introduction of the operator 0{A) makes it 
possible to identify a quantized action which is multiplicatively renormalizable, a feature 
which should not be related to a specific choice of the gauge fixing term Sgf, of course on 
the condition that the usual Yang- Mills action S, quantized using the gauge fixing Sgj, 
thus S = Sym + Sgf, is renormalizable. 

As we shall see, these conditions will lead us to consider the nonlocal gauge invariant operator 
of mass dimension two 

0{A) = -i / d^xF;, [{dYJ'f^^, . (1.3) 



Expression H1.3() can be made local by the introduction of a set of additional fields. Moreover, 
we will be able to prove that it is possible to identify a local and polynomial action which turns 
out to be multiplicatively renormalizable to all orders. 



The identification of this action and the algebraic proof of its renormalizability, as explicitly 
checked through the evaluation of the one loop anomalous dimensions, are the main results of 
the present investigation, signaling that the operator (|1.3j) could be relevant for a better under- 
standing of the issue of the gauge invariance of the dimension two gluon condensate. 
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Besides the renormahzability, we should also provide a suitable framework to discuss the possi- 
ble condensation of the operator (|1.3() . i.e. {0{A)) ^ 0, which would give rise to the dynamical 
gluon mass generation. Although being out of the aim of the present work, we remark that, in 
the Landau gauge, 5^^^ = 0, expression (|1.3|) reduces, to the first order, to the mass operator 



ah 



fiu 



I 



d'xA-A; 



+ higher order terms 



:l4) 



Thus, it is not inconceivable that a nonvanishing condensate (^A'^A'^j^ ^ might provide a sup- 
port in favor of a nonvanishing condensation of the operator ()1.3() . i.e. (F-X^F) ^ 0. 



The plan of the work is as follows. In section 2 we give an account of a set of nonlocal and 
gauge invariant mass operators which can be introduced in the Abelian case. These include the 
Abelian version of the operator 0{A) of eq. (|1.3|) . the operator A'^^^ recently discussed in OIEI) 
the Stueckelberg term as well as the nonlocal mass operator J d^xAj^Aj^^, where Aj^ stands for 

the transverse component of the gauge field Aj^ = ^6^i, — A^,. Interestingly, in the 

Abelian case, it turns out that all these gauge invariant operators can be proven to be classically 
equivalent, i.e. they reduce to the same expression when the classical equations of motion are 
used. In section 3 we present a detailed discussion of the non-Abelian generalization of these 
mass operators. We shall see that all operators introduced in the Abelian case possess a non- 
Abelian gauge invariant extension. However, the classical equivalence between them is now no 
longer valid. In particular, we point out that, in the non-Abelian case, the mass operator of 
ea. (|1.3() exhibits differences with respect to the operator A'^^^. As we shall see, the latter can 
be expressed as an infinite sum of nonlocal terms, a feature which makes almost hopeless the 
possibility of achieving a consistent localization procedure for a generic choice of the gauge fixing 
condition. Section 4 is devoted to the study of the localization procedure of the mass operator 
H1.3|) and of the rich symmetry content of the resulting action. In section 5, the identification 
of a suitable local and polynomial action is provided. Its multiplicative renormahzability in the 
class of covariant linear gauges will be established by means of the algebraic renormalization. 
Having developed the general properties of the mass operator, we devote section 6 to the com- 
putation of its anomalous dimension at one loop. Our conclusions are presented in section 7. 
For the benefit of the reader, we have found useful to collect in several Appendices the explicit 
derivation of some relevant features of the various mass operators considered in this work. 



2 Mass operators in the Abelian case. 

In this section we shall discuss a set of nonlocal gauge invariant mass operators which can be 
added to the Maxwell action 



1 



d^xF^^F^, . (2.1) 



Perhaps, the simplest way of introducing a gauge invariant mass term is through the nonlocal 
gauge invariant variable A^ 



Al = [6,. - ^ j • (2.2) 

Expression 1)2. 2|) is recognized to be the transverse component of the gauge field, df^Aj^ = 0, and 
is invariant under the gauge transformations, i.e. 

6Al = 0, (2.3) 
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with 

5A^ = -d^uj , (2.4) 
Thus, for the gauge invariant mass term one writes 



{A) = I d^xAlAl . (2.5) 



A second possibihty of introducing an invariant mass term is provided by the operator A^j^, 
which has been recently analysed in [SI IE1 • The operator ^^j^ is obtained by minimizing the 
quantity / d'^xA^A^^ with respect to the gauge transformations, namely 

02{A) = Ali^ = min j d'^xA^A^ . (2.6) 

Making use of the decomposition of the gauge field into transverse and longitudinal parts 

T _ (. 



it follows that 



I d^xA.A, = j d^xAlAl + I d^xA^^A^^ . (2. 



Observe that both terms in equation (|2.8|) are positive definite. Moreover, as discussed in 
the functional J d'^xA^^A^ achieves its minimum when d^A^j^ = 0, i.e. A^ = 0, so that 

Oi{A) = 02{A), (2.9) 

which establishes the equivalence between expressions (|2.5p and H2.6() . It is worth mentioning 
that the gauge invariant functional vl^j^^ has been proven to be an order parameter for the study 
of the phase transition of compact three-dimensional QED |S] . 

A third possibility of introducing an invariant mass operator in the Abelian case is by means of 
the Stueckelberg term |SH] 

03(A) = j d'^xiA^ + d^ct^f , (2.10) 

where <f) is a dimensionless scalar field. Expression (|2.1U|) is left invariant by the following 
transformations 

5A^ = -d^uj , 



The mass term (|2.10p can be rewritten in the form of a U{1) gauged fi-model, by introducing 
the variable 

U = e'^^ . (2.12) 
Thus ^ 

OsiA) = J d\ (^A^ - lu-%U^ . (2.13) 
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Transformations (|2.11j) read now 



U ^ UV , (2.14) 

with 

V = e'""^ . (2.15) 

One checks that the quantity — |C/~^9^C/) is left invariant by the transformations (|2.14jl . 
Analogously to the operator 02{A), expression H2.1U() can be proven to be classically equivalent 
to the mass term of equation (|2.5|) . This is easily seen by looking at the equations of motion 
which follow from the gauge invariant action 

S=\j d^xF^^F^, + 1^ j {A^ + d^(l)f , (2.16) 

namely 

d^F^^ + m^{A^ + d^(t)) = 0, 

d^(t> + d^A^ = 0. (2.17) 

In particular, from the second equation of H2.17() . we obtain 

<l) = -^dA, (2.18) 

so that 

O-M) = I d'x I^A, - ^aJ^ ' = I d'xAlAl . (2.19) 

Thus 

0^{A) = Oi{A), (2.20) 

which establishes the classical equivalence between expressions (|2.5() and H2.10() . Also, from 
H2.18() one sees that the scalar field 4> is related to the longitudinal mode of the gauge field j4^. 

Finally, a fourth mass operator can be introduced by considering the nonlocal quantity 

04{A) = -^j (fxF^^^F^, . (2.21) 

Again, this term is seen to be equivalent to expression (|2.5|) . In fact 

= ^jd^x A,^ {d^A, - d^d.A^) +A^^ {d^A^ - d.d^A, 

d'xA, [a, - ^a)j = j d^xAlAl , (2.22) 
thus 

Oa{A) = Oi{A). (2.23) 

Albeit nonlocal, the operator ()2.21|) can be made local through the introduction of suitable 
additional fields. More precisely, in the present case, one has 

-]rr? j d^xF^^^F^, ^ j d^x {^^,d^B^, + ^ {F^.B^, - F^.B^,)^ , (2.24) 
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where B^i, and B^i, are a pair of antisymmetric complex fields and m is a mass parameter. 
Eliminating Bf^i, and B^j^i, by means of their equations of motion, one gets back the nonlocal 
action (|2.21|) . One sees thus that, once cast in the local form, expression (|2.21j) looks renormal- 
izable by power counting. It turns out in fact that, in the Abelian case, the localized term in 
the right hand side of eq. (|2.24j) can be added to the usual QED Lagrangian without destroying 
its renormalizability. 

3 Mass operators in the non- Abelian case. 

As we have seen, there exist several ways of introducing nonlocal gauge invariant mass operators 
in the Abelian case. In particular, the four mass operators (|2.5|) . (|2.6|) . (|2.1fl|) and (|2.2H) turn 
out to be equivalent. Let us face now the more complex case of non- Abelian gauge theories. Let 
us start by considering the operator A'^-^^. 

3.1 Non- Abelian generalization of the operator A^^^^. 

The operator A'^-^^ of expression 1)2. 6|) can be generalized to the non- Abelian case by minimizing 
the functional Tr / d^x A^A"/^ along the gauge orbit of^^jMHEOlilllMlEHlElEl, namely 

= minTr / d^'xA^^A^^, 
m J 

A^ = v)A^u+-^v)d^u. (3.1) 

A few remarks are in order. Although the minimization procedure along the gauge orbit of A^ 
makes the operator ^^^^ gauge invariant, it should be underlined that the explicit determination 
of the absolute minimum achieved by the functional Tr j d^x A^A^ is a highly nontrivial task 
which, in practice, requires the resolution of the issue of the Gribov copies. It has been proven 
that the operator Tr J d'^x A^A^ achieves its absolute minimum along the gauge orbit of A^ 
[Sni IHUl 1^ inn inS]- Moreover, it is also known that, in general, it possesses many relative 
minima along a given gauge orbit. Therefore, one has to be sure that the correct minimum 
has been selected. This requires a detailed knowledge of the so called fundamental modular 
region, which is the set of all absolute minima in field space of the functional Tr J d'^x A^A^. 
The fundamental modular region is contained in the Gribov region, which is defined as the set 
of all relative minima of Tr J d'^x A^A^. While the Gribov region turns out to be still plagued 
by the presence of additional Gribov copies, the interior of the fundamental modular region is 
free from Gribov copies jSHl EOl 1^ 1121 IMl ISl lEj; a feature of primary importance for a correct 
quantization of Yang-Mills theories. However, a knowledge of the fundamental modular region 
of practical use in the Feynman path integral is not yet at our disposal. All this should give 
to the reader an idea of the real difficulty of obtaining an explicit expression for the absolute 
minimum configuration of the functional Tr J d'^x A^A^. A more modest program would be 
that of considering the Gribov region instead of the fundamental modular region, amounting to 
consider field configurations which are relative minima of Tr J d'^x A^A^. These configurations 
can be constructed in a relatively easy way as formal power series in the gauge field A^^. As 
discussed in Appendix a minimum configuration of Tr J d^x A^A^ is attained when u = h 
so that A|J is a transverse field, dfj_A^ = 0. The transversality condition can be solved order by 
order [H^, allowing us to express h as a formal power series in the gauge field A^, i.e. h = h{A). 
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This gives 



'flU 



Q2 



A. 



^9 



^9 
2 



(3.2) 



In particular, the configuration A^ turns out to be invariant under infinitesimal gauge transfor- 
mations order by order in the gauge coupling g |64j . see also Appendix namely 

SA^^ = 0, 
SA^ = 

Thus, from expression (|3.2|) it follows that 



-d^uj + ig [A^,i 



(3.3) 



Tr I d^A^A^^, 



^1 ( 



Q2 



a: 



gf 



abc 



dA'' 



Q2 



dA 



+ OiA") .(3.4) 



We see that the operator A^^^ can be expressed as an infinite sum of nonlocal terms. Such a 
nonlocal structure looks almost hopeless to be handled in a consistent way for a generic choice 
of the gauge fixing term. The only possibility here seems that of adopting the Landau gauge 
condition, d^A"^^ = 0. In this case, all nonlocal terms in the r.h.s. of equation (|3.4j) drop out, so 
that 

1 



4^ 



d X A. , , A. . , 



in the Landau gauge . 



It is worth remarking that, as proven in the massive Yang-Mills action 
, . . ...2 

5, 



\ j d'xF-,F-^ + ^ / d'xA^Al + 



d^x ( b^duAf, + ^daDf!'c^ 



(3.5) 



(3.6) 



where 6" is the Lagrange multiplier enforcing the Landau condition, d^A'j^ = 0, and c", c" are the 
Faddeev-Popov ghosts, is multiplicatively renormalizable to all orders of perturbation theory. 

In summary, we have seen that the operator ^^j^ can be generalized to the non-Abelian case. 
In addition, when treated as a formal power series in the gauge field A^, it has the pleasant 
property of reducing to the renormalizable operator J d^xA'^A'^ in the Landau gauge. 

We also recall that the operator J d'^xA'^A'^ turns out to be renormalizable to all orders of 
perturbation theory in the more general class of the linear covariant gauges [29j, a fact which 
has made possible to give evidence of a nonvanishing condensate (^A^A'^ in these gauges jl3| . 
However, outside of the Landau gauge, the relationship between A^j^ and / d^xA'^A'Jj^ is lost, 
so that a study of the nonlocal operator A^j^ becomes difficult. The operator A'^-^^ lacks thus 
a simple computational framework outside of the Landau gauge. 



3.2 Non-Abelian generalization of the operator J d'^xAj^Aj^. 

The discussion of the previous section allows us to generalize the operator f d'^xAj^Aj^ in the 
non-Abelian case. In fact, according to expression 1)3. 2|) see also Appendix^ it is possible 
to introduce a gauge invariant non-Abelian transverse field. It follows thus that the non-Abelian 
generalization of the mass operator f d'^xAj^Aj^ is provided by expression ()3.4() . This establishes 
the equivalence between the non-Abelian version of f d'^xAj^AJ^ and the functional A'^^^ within 
the space of the formal power series. Moreover, the operator f d'^xAj^Aj^ is plagued by the same 
difficulties affecting A^njjj. 
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3.3 Non-Abelian generalization of the Stueckelberg term. 

The Stueckelberg term, ea. (|2.1U|) . can be promoted to the non-Abehan case |58| . namely 

2 



with 



Os = Tt I (fx ( - -U-'d^U 



u 



(3.7) 



(3.8) 



where {T"^}, a = 1, ...,N'^ — 1, denote the hermitian generators of the gauge group SU{N), and 
where cp"" is a dimensionless scalar field in the adjoint representation. As shown in Appendix IbI 
expression (|3.7)1 is left invariant by the gauge transformations 



A 



uv 



The resulting non-Abelian massive action 

1 



Ss = -Tr / dt'xF^.F^, + m'Tr I d'^x { A 



(3.9) 



(3.10) 



looks local. However, it is not polynomial in the Stueckelberg field (j)"". In fact, when expanded 
in a power series in the field the term U~^d^U gives rise to an infinite number of vertices. 
This jeopardizes a consistent perturbative treatment of expression (|3.10() . To the best of our 
knowledge, the action (|3.1())1 is not multiplicatively renormalizable 58 , see also the recent 
discussion given in jSS]. As done in the Abelian case, it is interesting to have a look at the 
classical equations of motion which follow from the action (|3.1fl|) . i.e. 



D„ A„ - -U~^d,.U 



. 



(3.11) 



Equation (|3.1H1 can be used to express the Stueckelberg field 0" as a power series in the gauge 
field [nni, see also Appendix El yielding 



r'^^u^M^ + y^^'^dA'^dA' ] + o{A 



dA 



1 



pabc ' 



1 



Q2 



(3.12) 



Therefore 
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9 



r)A'^ 1 

+ ir''A%^ + -rbcdA'—dA' 



1 

Q2' 



+ 0{A^ 



(3.13) 



Thus 



Os 



Tr 


jd^x 




1 <fx 


2 J 




9f 


abc AaT 



1^ - -u-^d^u 



AfAf + 2gA 



dA"" 



1 



1 



dA^\ dA" 



f ^ [r'^Aid^^ + -r'^dA'-^dA" 

+ O(A') . 



(3.14) 
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Moreover, taking into account that, due to the transversahty of , the second term of the 
expression above vanishes by integration by parts, we obtain 



Os = ^Tr / d^x 



(3.15) 



which coincides precisely with expression ()3.4() . This shows the classical equivalence, within the 
space of the formal power series, between the Stueckelberg mass operator and the functional 



A^jj^ in the non-Abelian case. 



3.4 Non-Abelian generalization of the operator J d^xF^y-^F^y. 

It remains now to discuss the non-Abelian generalization of the operator J d'^xF^y-^F^y, a task 
easily achieved by replacing the ordinary derivative, d, by the covariant one, D, namely 



TV j d^xF^y^F,, = i 1 d^xF^y [{D-" 



-11 ^■^ 



(3.16) 



We remark that this term can be introduced in any gauge and, unlike the functional ^^in' 
does not require any specific knowledge of the properties of the Gribov region as well as of 
the fundamental modular region. It has already been considered in [HZI in the case of the 
three-dimensional Yang-Mills theories, where the use of the operator H3.16() was based on its 
appearance in e.g. the two-dimensional Schwinger model. However, so far, it has not yet been 
analysed in four dimensions. Although in the Abelian case the operator J d^xF^y-^F^y turns out 
to be equivalent to ^^in' ^^^^ more true in the non-Abelian case. This can be understood 
by observing that, thanks to gauge invariance, the expression ()3.4() for can be rewritten 

directly in terms of the field strength F^j^y. In fact, as proven in [^0], it turns out that 



>r / d^x Uy^^ 



■^F>i^Fi^y, J^DyFXfj^ 



-2^ ;^^am 



— n F \ —D F 



+ OiF^), 



(3.17) 



from which the difference between the operator H3.16() and Af^^^ becomes apparent. This in- 
teresting feature gives to the operator 1)3. 16(1 a privileged role with respect to the localization 
procedure. In fact, while in the case of A'^-^^ one has to deal with an infinite number of nonlocal 
terms, expression ()3.16() seems to be more manageable. In the next section the localization 
procedure of the operator (|3.16|) will be discussed. 



4 Localizing the mass operator Tr J d^xF^j^^F^^. 

The localization of the operator Tr J" d^xF^y^F^y can be achieved by generalizing the procedure 
adopted in the localization of the Abelian operator J d^xF^^y-^F^y, ea. ()2.24|) . Let us start by 
considering the Yang- Mills action with the addition of the mass operator (|3.16|1 . i.e. 



where 



Sym + So , 



Sym = ^ I d'^xF^yF^y , 



(4.1) 
(4.2) 
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and 



2 

m 



Sa = -'-fj d'xF^, . (4.3) 

The term 1)4. 3|) can be localized by means of the introduction of a pair of complex bosonic 
antisymmetric tensor fields in the adjoint representation, according to 

e-^o = J DBDB [detD^fe^^ [- Q j d^xB^^^OfD^fB^, + jd^^^^B- B)l^F-}j 

(4.4) 

where the determinant, (det D^) , takes into account the Jacobian arising from the integration 
over the bosonic complex fields (-B^i/, -B^;/) • This term can also be localized by means of suitable 
anticommuting antisymmetric tensor fields {G^^j^^G'^^) , namely 

{detD')' = I DGDGexp Q J d^xG^Df dI^G^^)j . (4.5) 

Therefore, we obtain a classical local action which reads 

Sym + Sbg + Sm , (4.6) 



where 



Sbg = \ j d\ (^B^^Dfot^B^^, - G^Dfo'fG';,,) 



ini 



Sm = — j d-x{B-B)l^F;,. (4.7) 

The localization procedure does not destroy the gauge invariance of the resulting action. In fact, 
it is easily checked that expression ()4.6() is left invariant by the gauge transformations 







b 




= gr'^u;' 


) DC 

B/^u 




= gr'^u;' 


B/^u 


^Gl, 


= gr'^u;' 




^Gl^u 


= gr'^uj' 





5(SyA/ + 5BG + 5™)=0, (4.9) 

so that condition (I.) is fulfilled. Let us proceed thus with the identification of a suitable quan- 
tized action, associated to expression (|4.6() . which enjoys the property of being multiplicatively 
renormalizable. For that, we follow the setup successfully introduced by Zwanziger j68U69j in the 
localization of the nonlocal horizon function implementing the restriction to the Gribov region 
in the Landau gauge. In a series of papers, Zwanziger has been able to show that the restriction 
to the Gribov region can be implemented by adding to the Yang-Mills action a nonlocal term, 
known as the horizon function, which is given by 

^Horiz = l^g" I d^xr^'Al [M-^y f^'Al , (4.10) 

where 7 denotes the Gribov parameter |70| and Jv[°-^ is the Faddeev-Popov operator of the 
Landau gauge 

A^"" = -d^ [d^S'^' + gr'A^) . (4.11) 
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As proven in j68[ I69j . the nonlocal horizon term (|4.1()j) can be localized by means of a suitable 
set of additional fields, in a way analogous to that of ea. ()4.4l) . Remarkably, the resulting theory 
is renormalizable to all orders, obeying the renormalization group equations. Thus, it seems 
natural to us to adopt here the same procedure. According to [^11101]) we treat the operators 
^Imu^^lu ^liv^llvi entering the expression for Sm in eq. (|4.7j) . as composite operators coupled 
to suitable external sources V"o-p^,y(a;), Vap^v{x). This amounts to replace the term Sm by 

\ I d'^x [V„p,,B-^F-^, - V^p^^B^^^F-,) . (4.12) 

At the end, the sources ^apfiui^^i ^appui^x) are required to attain their physical value, namely 



^appv 



— Vcrppu 

phys 



— ITTl 

= -TT {^crp^pu - 5av5pp) , (4.13) 
phys Z 



SO that expression ()4.12|) gives back the term Sm- As pointed out in |SSlinn]) this procedure allows 
us to study the renormalization properties of the Green's functions obtained from the action 
{Sym + Sbg) with the insertion of the composite operators B^^^F^^ and B'j^i^F^^^. Following 
[HHlinni) let us focus first on the properties of the action {Sym + Sbg) which, as we shall see, 
displays a rich symmetry content. 

4.1 BRST invariance. 

In this section we shall discuss the symmetry content of the action {Sym + Sbg)^ where Sym is 
the Yang- Mills action, eq. (|4.2|) . and Sbg depends on the localizing fields (S^^, 5^^,, G^^, G^j^) , 
ea. ()4.7() . Let us begin by introducing the gauge fixing term, chosen here to be that of the linear 
covariant gauges, namely 

S = SYM + SBG + Sgf , (4.14) 

with 

5,/ = I <fx + h'^dpAl + c"5^<c^) , (4.15) 

where 6*^ is the Lagrange multiplier and c", c" stand for the Faddeev-Popov ghosts. It turns out 
that the action S is left invariant by the following BRST transformation, i.e. 









sc" 


= Ir'^c-c^ 




^^p,i/ 


= gr'"^c'B-^, 






= gr'^c^B^, 


; 






) 








S(f 


= h\ 




sb" 


= 0, 




^2 


= 0, 





(4.16) 
and 

s5 = . (4.17) 

This is easily verified by observing that the term Sbg can be written as a pure BRST variation, 
according to 



1 



Sbg = -^s I d''xG1,DfD''^-B-^^ . (4.18) 
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fields 


A 


c 


c 


h 


B 


B 


G 


G 


dimension 


1 





2 


2 


1 


1 


1 


1 


ghost number 





1 


-1 











1 


-1 


Q /-charge 














1 


-1 


1 


-1 



Table 1: Dimension, ghost number and Q/-charge of the fields. 



Of course, the same property holds for the gauge fixing term 



5, 



9/ 



: j d^x i^^h'' + (fd^Al^ . (4.19) 



In addition to the BRST invariance, and in complete analogy with the Zwanziger action |681l69j 
implementing the restriction to the Gribov horizon, the model displays a global invariance U{f), 
/ = 6, expressed by 

Qi..a^S = , (4.20) 

where 

The presence of the global invariance U{f) , / = 6, means that one can make use of the composite 
index i = {/xi^}, i = (1, . . . , 6). Therefore, setting 



f T3<^ r'"^\ ^ f R'* \ 00^ 



we get 



Sbg = I d'x (BfDfDl^Bt - G^DfD';Gt) , (4.23) 
and for the symmetry generator 

By means of the diagonal operator Q/ = Qu, the i— valued fields turn out to possess an additional 
quantum number, displayed in Table Q together with the dimension and the ghost number. 
Besides the global U{f), / = 6, invariance, the action (|4.14j) possesses the following additional 
rigid symmetries 

7^Jf 5 = , (4.25) 

where A = {1, 2, 3, 4} and 

= Id'xi BfA: -G'i ^ ^ 



J \ ' 5G": ^ 5Bf 
nf} = I d'a:iB?S, + G^ ^ 



■i] \ ' 5G": ^ 5Bf 



(4.26) 
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sources 


U 


V 


u 


V 


dimension 


1 


1 


1 


1 


ghost number 


-1 





1 





Q /-charge 


-1 


1 


1 


-1 



Table 2: Dimension, ghost number and Q/-charge of the sources. 



Let us conclude this section by showing that also the source term (|4.12|) can be introduced in a 
BRST invariant way. This is achieved by considering the following source term 

+ X2Ui^,yd^daViua — C [Ui^iyVif^i^VjapyjalS — Ui^^Vi^i^UjapU , (4.28) 

with 









= 0, 








= 0, 




= . 



(4.29) 

The quantum numbers of the sources are displayed in Table [21 Therefore, for Saux one gets 
Saux — d X \Ui^iiG ^ifiuB^ Vi^i/B^ -\- UijipG^ F^y -\- xi (y^i^ud ^i^u 

+ Vi^iuVifiuVjalSyjaP — '^Ui^uUi^uVjapVjap)'] ■ (4.30) 

The parameters xi; X2 and C ai"6 free parameters, needed for renormalizability purposes. The 
action Saux reduces to the term Sm of ea. H4.7|) when the sources (Vifj,y,\^iixuiUi^y^Ui^y^ attain 
their physical values, given now by 

{y^illV t^illV ^Ui^y ^Ui^y^ — (Va-p^y ,Vcrp^y ^ Uo-p^y , Ufjp^y^ , (4.3l) 



phys 



phys 

Uap^y = U^ppu = 0. (4.32) 



Thus 



S, 



aux 

phys 



Sm-^ I d^x Cm^ , (4.33) 



so that the term Sm is recovered, modulo the constant quantity Qrn^. All ingredients needed to 
study the renormalizability of the action 

Sym + Sbg + Sgf + Saux , (4.34) 
are now at our disposal. This will be the task of the next section. 
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sources 


n 


L 


Y 


Y 


X 


X 


dimension 


3 


4 


3 


3 


3 


3 


ghost number 


-1 


-2 


-1 


-1 


-2 





Q /-charge 








-1 


1 


-1 


1 



Table 3: Dimension, fermion number and Q/-charge of the external sources. 

5 Identification of a multiplicatively renormalizable action. 

In order to discuss the renormalizability properties of our model, we have first to write down 
all possible Ward identities expressing the symmetry content of the starting classical action, 
eq. (|4.34)) . Let us begin by working out the Slavnov- Taylor identity. Following the algebraic 
renormalization procedure as described in [T^, we need to introduce additional external sources 
(OJ^, L"", Y-^, Y^^Xf, Xf^ in order to define at the quantum level the composite operators entering 
the nonlinear BRST transformations of the fields {A^^, , , , , G'^) , eqs.^JE^. In the 
present case, this term reads 



with 
and 



n'^^A'', + L'^c" - 


- Y^-Bt - YtBt + X'^G'i + XtG'i) , 


(5.1) 




sL'^ = , 


(5.2) 


sYt 


= 




sXf 


= 


(5.3) 


sXt 


i 


(5.4) 


sYt 


= . 


(5.5) 



The quantum numbers of the external sources L°, 1^", 1^", , X?) are displayed in Table 
131 For the complete action S 

S = SyM + Sgf + SbG + Saux + Sext , (5.6) 

we obtain 

S = Sym + j d*x (^^h%'' + h''d^A1 + <fd^D^^c^'^ + j d''x{B'^DfDl'B',-G'iDfDl'Gi) 
+ j d^x {{Ui^^G^ + Vi^,Bt - Vi^,Bf + U^^^Gt) F^, + xi {%ud^Vi^u - %ud^Ui^,) ) 

+ J d'^XX2 {yifj.udfj,daViua — Ui^ud^daUiua) — J d'^x( [Ui^^Ui^^UjapUjalB + yifiuVifiuyjalSyjalS 

- 2%,U^^,V,^pyj^fs) + I d^x (-%Dfc'> + ^-r'^L'^c'c'^ + gr'%-c'>Bt + gr'^Y^c^'Bt 
+ gr'^Xtc^Gi + gr'-Xfc'ci) . (5.7) 
Expression (|5.7j) obeys several Ward identities, which we enlist below 
• the Slavnov- Taylor identity 



5(S) = 



(5.8) 
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+ 



(5S 



5Bf SY'^dBf 5Xf5Gf 



5f 1 ^ + Vi^^-^ + UiuuTT-. Y,^-^ + Xf- 



6X? ' J 6Gf 
• the global U{f) invariance , / = 6, i.e. 



(5.9) 



QijS = 



(5.10) 



where 



u. 



+ V, 



SX" 



6Xf 



the exact rigid symmetries 



where A = {1, 2, 3, 4} and 



7^: 



(4) 



d^x Bf 



d^x Bf 



5G^ 



5G'^ 



5G'^ 



d^x Bf- 



5G^ 



n[fj: = 



(5.12) 



U- ^ I y" I .Y" ^ 



G'iJ^ + Vi„,,^J— + U, 



' 6Xf ^ ^ 6Yf 



r^'^-v- ^ I u — 



G 



^ 5Bf 



Vi 



6 



'^"5Vi^.. ' bX-^ nY- 



(5.13) 



the gauge fixing condition 



the antighost Ward identity 



where 



= 



(5.14) 

(5.15) 
(5.16) 



5.1 Determination of the most general local invariant counterterm 



Having established all the Ward identities fulfilled by the complete action S, we can now turn 
to the characterization of the most general allowed counterterm Yf. Following the algebraic 
renormalization procedure |71j . is an integrated local polynomial in the fields and sources 
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with dimension bounded by four, with vanishing ghost number and Q/-charge, obeying the 
following constraints 



ST." 



= 0, 
= 0, 

= 0, 
= 0, 



in addition to 



Bj:^^ = , 



(5.17) 
(5.18) 



where Bs is the nilpotent linearized Slavnov- Taylor operator 



0s 



5Y S (5E (5 5Y 5 (5E (5 S 



ST. 6 5T S 



f 5T 



+ 



5Bf SYI" 5Y'' 5Bf \6Bf 



5 5T 5 



(5y" 5Xf5Gf 



SGf ' / 6X? 



6a 



HIV _ 



Bj^Bj: = . 

After a rather lengthy analysis, for the most general allowed counterterm we have found 



(5.19) 
(5.20) 



E*^ = aoSvM + ai d xA 



5Sy 



M 



+ j d^x ((ai + a2) (0» + 5^cf ) d^d' + as^/"'^ K + 5^^) - a2p''''L-c'c') 
+ J d^x{ (2a3 + 04) B^d^Bt - (2a3 + 04) G^^^Gf 

- (ai + 2a3 + 04) gf'Bf {d^A'^ + 2A%) Bf + (2ai + 2a3 + 04) f"'' f''' Bf Af^A^ 



+ (ai + 2a3 + 04) gr^Gf [d^A;, + 2A'^d^) G\ - (2ai + 203 + 04) 5^"" V'Gf^^^G; 
- a^gr^'c^ {y^Bt + Y^Bt - X\G\ - X\Gi) + [(ai + as + as) Id^A^ 
+ (2ai + 03 + as) (t^^M-Gf + Vi^^B^ + C/,^,Gf - %^.B^) 

yahcd 



+ 
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~ GIG, I j [BjBj — G^jGj \ + a-j {B^B^ - G^G^) {Vi^^Vi^v — Ui^yUi^v) 



+ 0,8 i^B"^ G'jVifj^yU j + Gf G'jUifj^yUj fj_i, + B^ BjVif^yVj ^y — GlB'jU-myVjpi,y 

1__ 1-- \ - - 2 

~ 2^i^'j^it^vyjnu + —G'^GjUifj^yUj^u j + agC {Vi^yVi^y — Ui^j^yUi^j^J) 

+ oioXi (Vi^i/SVi^,, - t/j^j.5^L/i^,,) + aiixi (t^i/.j/5M^a^iiva - ^^^^ , (5.21) 

where (ao, ai, a2, 03, 04, 05, ay, ag, ag, aio, an) are free parameters and X"-^'^'^ is an invariant tensor 
of rank four with indices in the adjoint representation and such that 

^abcd ^cdah 

^abcd ^ ^bacd _ (^5 22) 
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For a general discussion of the properties of higher rank invariant tensors we refer the reader to 
|72j . Let us only mention that an invariant rank 4 tensor like A"^'^'^ obeys a generalized Jacobi 
identity 

jman ^mbcd _j_ jmbn ^^amcd _|_ jmcn yabmd _j_ jmdn ^^abcm q 

These parameters aj,i = 0,...,ll, should correspond to a multiplicative renormalization of the 
fields, parameters and sources of the starting classical action S. However, it turns out that the 
counterterm (|5.21|) cannot be reabsorbed through a renormalization of the parameters and fields 
of S. This means that the starting action S is not stable against radiative corrections. Said 
otherwise, S is not the most general local invariant action compatible with the Ward identities 
(|5.8|) - (|5.15|) . In fact, from the expression (|5.21|) it follows that the term 

/r \abcd , _ _ \ / _ _ 

d^x a, - GfGt) {V,,,V,^, - U,,,U,,,) + [sfBl - G'iG\) [B^Bf - G^j 

+ ag [B'^G'jVi^uUj fj,u + G^G'jUifj^uUj^u + B^ B'^Vi^yVj ^i, — G^B'-Ui^yVj^u 

— G1 B'jUifj^uVj fjLu + G1 B'jUifj_i,Vj — —B^BjVifj^i^Vj^i^ + -^GIG'jUi^uUj^v 



~ 2^iBjVi^uVj^i, + -GiGjUi^uUjfj,J^ 



(5.24) 



fulfills all Ward identities. Moreover, this term does not correspond to a renormalization of the 
parameters and fields of S. This follows by noting that the counterterm (|5.24j) is in fact absent 
in the expression ()5.7|) . 

A stable action S is thus obtained by adding to the action S the following expression 
5a = y d^x [Ai {BtB^ - GtG'l) {Vj^.Vj^, - U^^^Uj^,) 

\abcd , _ _ \ / _ _ 

+ ^ ( BfB^ - GfG^ ] ( WBf - G'iG' 
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J J 33 



+ A3 [B'j^GjVi^iyUjfj^i/ + GlG'jUi^yUj^y + B^ B'jVi^yVjp^y — G^ BjUi^iyVjfj^i/ 

-B°-BjVifj,uVj^u + -< 



Bj Uifxi/Vj ^1/ ~\~ Gj Bj UifiiA/^j fjLU „ ^ifMU^j ^1/ ~l~ „ Gi Gj UifjiuU j 



1 -a -a 1 -a "a 

-BfB'-Vi^yVj^y + -G°;G'jUifj,uUj^u 



(5.25) 



where Ai,A3, are free parameters, namely, by taking as starting point the action 

S = SyM + Sgf + SbG + Saux + S\ + Sext ■ (5.26) 

The previous algebraic analysis can be repeated for the action S. For the most general allowed 
counterterm we find now 



18 



= aoSvM + a^ j d^xAl^-j^ 

+ j d^x ((ai + a2) + 5^c«) a^c'^ + 025/'^"'= + d^<f) - a2p''''L-c'c'^) 

+ j (fx{ (2a3 + 04) B'id'^B'i - (2o3 + 04) Gf^^G" 

- (ai + 2a3 + 04) ^/"''^A" + 249^) + (2ai + 2a3 + 04) 5 V'/^'^'A^^^^^i 
+ (ai + 2a3 + a4) ff/^'^Gf (^^4 + 2A^a^) Gl - (2ai + 2a3 + 04) ^'/"'''/'''GfAjA^G,^ 

- a25/"'^c« fF/Sf + y/Bf - Xi-Gf - X^Qi] + [(ai + as + as) 29^^" 



+ (2ai + as + as) ^/"''^^^^J (t/i^-Gf + F^^^Sf + Ui^.G'i - Vi^^Bf) 

\abcd . _ _ \/_ _ \ 

+ (4as + Se) ^ [BtB\ - G1G\) [B^Bf - Gp'^) 
+ (2as + 57) Ai [BfB^ - G^G^) (F.^.y.^. - a,^,U,^,) 

+ (2as + as) A3 {BfGjVi^^Ujfj,^ + G^GjUifj^^Uj^i, + BfBjVifj_,^Vjfj_^ — G°5j?7i^i,1^-^i, 

1__ 1-- \ - - 2 

~ ^^i^'j^ilJ-v^m^ + 2^iG'jUi^uUjni, 1 + agC {ViixuVifxij — Uini/Uinv) 

+ aioXi - Ui^yd'^Uii^y) + aiixi (Vi/i!/5^«9Q^i;/a - Uifj,udfj,daUii,a)} , (5.27) 

As a useful check, let us show that S*^ can be reabsorbed by means of a multiplicative renormal- 
ization of the parameters, fields and sources of S. Setting 

Jo = ZjJ , 

6 = Z^^, (5.28) 



where 



(p = {A,b,c,c,B,B,G,G} , 

J = {n,L,U,U,V,V,X,X,Y,Y} , 

e = {<?,«, XI, X2,C,Ai,A«''^^ As}, (5.29) 



it follows 



S((/)o, Jo, Co) = ^(P, J, + V^'{<P, J, + 0{v^) . (5.30) 
In particular, the renormalization constants are found to be 

Zy^ = 1-^(1 + 1) , (5.32) 

Zl/' = Zy^ (5.33) 

^1/2 ^ ^-1/2^ ^534^ 
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= 


(5.35) 


Zl 


= Za , 


(5.36) 


^1/2 

Zb 


^1/2 ^1/2 rzl/2 , , / 

= =^G = = ^ + ^ ("3 + y j , 


(5.37) 


Zv 


= Zy = Zu = Zu = l-r](^ + ^ -a5) , 


(5.38) 


Zv 


X Y c ^A B ' 


(5.39) 


Zg 


2 ' 


/ r Ar\\ 

(5.4Uj 


Za 




(5.41) 


Zx, 


= 1 + e (ao - 2a5 + 07) , 


(5.42) 


^abcd 


= 1 - e (204 - ae) , 


(5.43) 


Zxs 


= 1 + e (ao - 2a5 + as) , 


(5.44) 


Zxi 


= 1 + e (ao + a4 — 2a5 + aio) , 


(5.45) 


Zx2 


= 1 + e (ao + a4 - 2a5 + ail) , 


(5.46) 


Z( 


= 1 + e (2ao + 204 — 4a5 — ag) . 


(5.47) 



5.2 Summary. 

In summary, we have been able to identify a local and polynomial action, given in expression 
(|5.26|) . which displays multiplicative renormalizability. This has been achieved by adding to the 
action S the term 5";^, ea. (|5.25|) . which is compatible with the complete set of Ward identities. 
When the sources {Vifj_u,Vifj,u,Ui^u,Ui^u) attain their physical value, ea. (|4.32() . Sx becomes 



Sx 








phys 





-m^Ai {B'^^B^^ - G^j^G^j,) + 



m 



2^ 
32 



+ 



1 abed 



16 



73a Tjb 



yoa /^b \ I rye jjd f^c /^d 



(5.48) 



This expression reminds us of a kind of Higgs term. There are, however, several differences. 
These are due to the antisymmetric character of the fields (-B^;^, S^^, G^j^, G^^,) with respect 
to the Lorentz indices. Moreover, we remark that, while {B'^^, B'^^^ are bosonic, the fields 
(G^^,G^j,) are anticommuting. With the exception of the term containing the parameter A3, 
expression (|5.48() displays thus a supersymmetric structure, a feature supported by the fact that, 
according to (|4.16|1 . the fields (-B^i/, -B^jy, G^^, GJJ^) transform as BRST doublets. Therefore, 
a certain number of cancellations among the contributions arising from these fields might be 
expected in the evaluation of the Green's functions of the model. The possible use of this su- 
persymmetric structure will be explored in the future, as well as its possible consequences for 
the Green's functions of the model. 



To conclude, let us give explicitly the starting action when the sources {Vif^y, Vifj,^, Ui^j,y, tJi^j,^) at- 
tain their physical value, eq. (|4.32|) . while the additional external sources (17^, L", , ,Xf,Xf^ 
are put equal to zero. 
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Sym + Sbg + Sm + Sx 
1 im 



phys 



1 



2 



+ 



I ahcd 



16 



fiiy^liiy ^fii/^fiuj \^f)a^pa ^ pa^ per 



m..BL - GLG?, J (BLBt - GLGl) + ^b%'' + b^d^Al + c^S^Z?!^''-^ 



■'p^p 



(5.49) 



Let us finally notice that each of the terms in ea. H5.48p is invariant w.r.t. to the gauge transfor- 
mations 1)4. 8() . More precisely, one has 



^ Sym + Sbg + Sm + S^ 



(5.50) 



phys , 



6 One loop renormalization. 



We now turn to the details of the explicit one loop renormalization of the Lagrangian (|5.49|) in 
the presence of the nonlocal operator. It is first worth noting some of the key features of ()5.49p 
in relation to the extraction of the one loop renormalization constants prior to discussing their 
calculation. First, considering the case when m is zero then one has a gauge theory fixed in an 
arbitrary linear covariant gauge where in addition to the usual gluon and Faddeev-Popov ghost 
fields there are two additional auxiliary fields, B^^, and G^j^ where the latter is anticommuting. 
Since these fields originate in localizing the nonlocal operator, when that operator is absent at 
m = 0, these new fields ought to play a completely passive role in the (one loop) renormalization. 
In other words the gluon, Faddeev-Popov ghost and quark renormalization constants ought to 
be equivalent to those obtained when and GJ^^, are formally absent. However, when they 
are present the algebraic renormalization formalism has demonstrated that they generate a new 
quartic interaction through (one loop) renormalization effects* which is indicated by the term 
with the independent coupling X"-^'^'^ in 1)5. 49() . In other words if one computes the B°^B^B'^B'^ 
four-point function at one loop with A"^'^'^ initially zero, there will be a divergent contribution 
at 0{g^) which will be removed by the counterterm generated by the term involving A"'"^'^. This 
is akin to the situation in A(/i>^ theory where the Lagrangian is multiplicatively renormalizable 
in four dimensions. However, the interaction can be replaced by a cubic vertex involving an 
auxiliary scalar field. The renormalization of this version of the Lagrangian still proceeds as 
usual except that the Lagrangian ceases to be multiplicatively renormalizable since a (j)'^ vertex 
will naturally be generated from one loop box diagrams. The standard A(/>^ /3-function and 
renormalization group functions can still be extracted with the auxiliary field version but one 
has to take account of the effects of the generation of the extra interaction. Indeed a similar 
situation arises in two dimensional four-Fermi theories where a formalism has been developed 
j73j and used to perform three loop calculations. The situation for our current Lagrangian is 
the same. The quartic interaction is generated via loop interactions and will be 0{g^). Thus it 
does not need to be taken into account for the extraction of the one loop anomalous dimensions 

*It might be useful to remark here that, at one loop order, the invariant rank four tensor which emerges from 

explicit calculations turns out to be proportional to g'^ ^f^av jebq ^mcv ^mdq ^ J^eap j^ebq j^mdp j^mcq^ ^ ^y^^^y^ j^jgjjg 

fact the conditions 15.221 1. 
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we are interested in. 



For the case when m is non-zero, there is a similar situation. The algebraic renormalization 
demonstrates that the now localized mass operator {B'^^, — B'^j^)F^j^, which is dimension three, 
mixes into two gauge invariant dimension two operators being those associated with the cou- 
plings Ai and A3. In other words computing the renormalization of the operator with a massive 
gluon propagator will inevitably lead to the generation of these two additional operators. As 
such this is nothing new in that it follows the pattern already known for the renormalization 
of local composite operators (See, for example, [H]). Indeed it is reassuring that this property 
emerges in an elegant way from the algebraic renormalization formalism for a localized nonlocal 
operator. However, these two additional operators do not form the complete basis of the possi- 
ble dimension two operators that higher dimensional operators can mix into when one uses the 
massive theory. Since each combination of pairs of the set {B^^^, B'^^, G^^,, GJJ^} are individually 
gauge invariant operators, to correctly treat the renormalization one would have to construct 
the full mixing matrix for this set. Though only those combinations with zero ghost number 
would be of importance. As we are primarily focused on extracting the anomalous dimension 
of the nonlocal operator itself, it will be apparent that this mixing matrix is not immediately 
required and we will defer its computation to a later article. 

Having outlined the status of (|5.49j) it is now evident how one goes about extracting the renor- 
malization constants which will lead to the anomalous dimension of F°:,,^F^„. Since we have 
localized this operator to {B'^^ — B^j^)FI^^ then the anomalous dimension of F^^^^F^,^ is equiv- 
alent^ to that of the gauge invariant operator {B^j^^ — B^j^)F^jj. Therefore, we can extract the 
anomalous dimension by inserting (-B^j/ — -B^i/)-?"^!/ into a -B^^vl^ two-point function and compute 
it using massless propagators. This is similar to how one determines the quark mass anomalous 
dimension by inserting the mass operator ipip into a quark two-point function, j75| ITfij . For 
— we will need the -6°^, anomalous dimension. However, we have carried out the 

full renormalization of all the fields of (|5.49|) at one loop by making use of symbolic manipu- 
lation programmes. The Feynman diagrams for the relevant Green's functions are generated 
with the Qgraf package, [ZZI) converted into Form, |^, input notation before extracting the 
divergences with the Mincer package, j79j . This uses dimensional regularization in d = 4 — 2e 
dimensions and we will remove the infinities with the (mass independent) MS renormalization 
scheme. If we define 

d 

70(a) = fi—lnZ^, (6.1) 

for (j) € {A1j^,c'^,ip,B'^j^,G'^^j} where a = (7^/(167r^), then the renormalization constants give the 
explicit results 

7^ (a) = [{3a-U)CA + 8TFNf]^ + 0{a^) , 

7e(a) = (a-3)CA^ + 0{a^) , 
7^(a) = aCpa + O(a^) , 

7ij(a) = 7G(a) = (a - 3)C^a + 0(0^) , (6.2) 

where Nf is the number of quark flavours^ , T'^T'^' = Cpl, P'^'^f"'^ = Ca^'"'' and tr (T^T^) = TpS^'K 
For completeness we note that the massless momentum space propagators of the fields are 



^Up to an overall scaling factor. 

■'■Although we did not consider matter fields in the previous analysis, it turns out that the multiplicative 
renormalizability of the action E, eg. 115. 26^ 1. can be extended to the case in which spinor fields are present. 
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-P)) = 


(c"(p)c-^(- 


-P)) = 


(V'(p)^(- 


-P)) = 




-P)) = 




-P)) = 



p2 



[I- a 



,PfiPu 



gab 



gab 

2p' 



(6.3) 



where p is the momentum. It is worth noting that the expressions for the gluon, Faddeev- Popov 
ghost and quark are equivalent to those obtained in the absence of and G^^, as expected. 
Indeed from examining the contributions from the diagrams involving these fields it is evident 
that the anticommuting property of introduces the necessary minus sign to exactly cancel 
the contribution from the graph involving a B'^^, loop. To verify that the and anomalous 
dimensions are correct, aside from correctly satisfying the equality demanded from the algebraic 
renormalization, eg . (15.37)1 . we have also renormalized both the gluon- and gluon-G vertices at 
one loop and verified that the correct one loop a independent coupling constant renormalization 
emerges as 

1^ + 0{a^) . (6.4) 



/3(a) 



Hence the renormalization of the operator -B^j^F^j^ proceeds by inserting B^^F^^ into a gluon- 
two-point function and extracting the divergence from the five one loop diagrams. Although we 
are regarding B^^^FI^j^ as multiplicatively renormalizable, since it is of dimension three it could 
in principle mix into the dimension three quark mass operator, ipip. However, at one loop there 
are no mixed diagrams of inserting ipip into a gluon-i? Green's function or of inserting B^ 
into a quark two-point function. If we define 



pa 



Oo = ZoO, 
where the subscript o denotes the bare object, with 



O 



r>a pa 



then we find 



Hence, with 



we deduce 



Za 



1 + 



11 



-Ga 



+ 0(a^ 



d 



70 (a) = fj,— In Zo 



70 (a) 



a + 0{a^) . 



(6.5) 

(6.6) 
(6.7) 
(6.8) 
(6.9) 



As the original operator was gauge invariant it is reassuring to note that 70(0) is independent 
of a. It is worth underlining here that the anomalous dimension 70(0) is equivalent to the one 
loop /3-function, where the overall factor of 2 is accounted for by noting that this is equivalent 
to the anomalous dimension of m as opposed to that of m^. This is interesting for various 



23 



reasons. First in the one loop renormalization of two-leg higher dimension operators in Yang- 
Mills theories the operators F^j^F^^^, D^F^^D^F^^ and D^D^F^pD^DyF^p each have the same 
one loop anomalous dimensions which is also the /3-function^ , |8H I82j . What is intriguing 
in the present situation is that the nonlocal operator F^^^^F^^^, which has a similar Lorentz 
contraction as the higher dimension operators noted above, has an anomalous dimension which 
is the same at one loop. There would appear to be no a priori reason either from the algebraic 
renormalization or other methods to expect this. Obviously, having the two loop correction 
to (|6.9|) would enhance our understanding of both the renormalization and significance of this 
nonlocal operator. With the exception of F^^^F^j^, the renormalization group behaviour of the 
two-leg higher dimension operators is also unknown^. It would be interesting to pursue this 
study to find out if a gauge invariant and renormalization group invariant mass dimension two 
condensate could be found using F^^^^^F^^, provided the operator condenses. Since evidence 
for the existence of a non-zero dimension two condensate arises in the fitting of data for gauge 
variant objects ^112011^1231151) as a first step, it would seem natural in the light of (|6.9|) to 
find out whether one could extract an estimate for the one loop renormalization group invariant 
condensate {dsF^y^F^j^) by fitting for 1/Q^ power corrections in measurements of correlations 
of gauge invariant operators. We refer to [841 1851 186j for a review of the role of such 
corrections which go beyond the standard SVZ-expansion |571 |HH1 EH ■ 



7 Conclusions. 

In this work the properties of the nonlocal gauge invariant operator Tr j d'^xF^y{D'^)^^ F^y of 
mass dimension 2 have been investigated. We started by looking at the Abelian case, where 
several nonlocal gauge invariant operators have been considered. Moreover, in this case, all 
operators turn out to reduce to the same expression when the classical equations of motion are 
employed. All Abelian operators generalize to the non-Abelian case. However, their classical 
equivalence does not hold anymore. In particular, the operator Tr J d'^xF^y{D'^)^'^F^y exhibits 
differences with respect to the operator A^^^. 

Albeit nonlocal, the operator Tr j d'^xF^^^^D^)^^ F^^, can be cast in local form by the intro- 
duction of a suitable set of additional fields, in contrast with the operator j4^jj^. A local and 
polynomial action has been identified, eq. (|5.2())) . and proven to be multiplicatively renormalizable 
to all orders in the class of linear covariant gauges by means of the algebraic renormalization. We 
point out that this action possesses a finite and relatively small number of parameters, a feature 
useful for higher order computations. We have calculated the one- loop renormalization group 
functions of the model. We have recovered the anomalous dimensions of the elementary fields, 
if already known. In the case of the nonlocal operator, we have found that the renormalization 
group behaviour is dictated by the /3-function at one- loop. 

The possibility of having at our disposal a local and renormalizable action might provide us 
with a consistent framework for a future investigation of the possible existence of the conden- 
sate (F-^F). 

^We have checked that the dimension ten operator D^D^DaFpgDfj^D^DaFpg has the same one loop anomalous 
dimension too, 80 . 

^For details concerning the renormalization (group) properties of F^^F^j, with or without massless/massive 
quaxks, we refer to |83|. 
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A Properties of the functional /a[w]. 

In this Appendix we recall some useful properties of the functional fA[u] 

fA[u] =Tr j (fxAlAl = Tt j (fx [v) A^u + -^v)d^v^ {u^ A^u + ^u^d^u^ . (A.l) 

For a given gauge field configuration A^, /aIu] is a functional defined on the gauge orbit of A^. 
Let A be the space of connections A^ with finite Hilbert norm ||j4||, i.e. 

1 1^1 1' = Tr y d^x A^A^ =^ J d^xA^A^ < +00 , (A.2) 

and let U be the space of local gauge transformations u such that the Hilbert norm ||u^9u|| is 
finite too, namely 

lln^a-ulp = Tt J d'^x (u'^d^.uj (^-u^a^n) < +00 . (A.3) 
The following proposition holds ISHl IS3 EH E21 E31 
• Proposition 

The functional /aIu] achieves its absolute minimum on the gauge orbit of A^. 

This proposition means that there exists a h ^li such that 

5fA[h] = 0, (A.4) 
&^fA[h] > 0, (A.5) 
fA[h] < /aM, yueU. (A.6) 



The operator A^-^^ is thus given by 

Al,, = minTr / d^x A^^A; = /^[/i] . (A.7) 

Let us give a look at the two conditions HA.4|) and HA.5|) . To evaluate 5fA[h] and 5'^fA[h] we 
set" 

V = he^3^ = he'S'^"^^ , (A.8) 



rpa rpb 



" The case of the gauge group 517 (A'') is considered here. 
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where oj is an infinitesimal hermitian matrix and we compute the linear and quadratic terms of 
the expansion of the functional /^b] in power series of u. Let us first obtain an expression for 



A'" 



Al = v^Ai^v + -v^df,v 
M ^ g ^ 

g V ^ y g 

'^9 



(A.10) 



Expanding up to the order cj^, we get 



A'" 



igu- g 



A^^ (1 + igu - g'^%-] +-{l-iguj- g'^%-] 9. (1 + igu: - g 



(l - igu - (4 + - 5^4^^) + 



+ -(1-iguj-g— ][ igd^^uj - — (5^u;) u - —uj {d^ 



UJ 



4 + i^^o; - ^Ay - zguAl + g^uA^u - ^-u^A^ 



+ - igd^j^u) - — {d^u) u - —ujd^uj + g ud^u) + 0{u-^) 



(A.ll) 



from which it follows 



Al = Al + ig[AlM + U[u,A%u]-d^u + i%,d^u]+0{u^) , (A.12) 



We now evaluate 



Tr J d^xA'^.A 



= Tr Jd^x (^Ayig[A%u] + ^[[u,A%u]-d^u + i^[u,d^u] + 0{u' 



A';, + ig[A';„ CO] + ^[[u, a;:], u;] - d^u + i^[u, d^u] + 0(0 

+ i^A^, d^u] + ig[A% - g^ f^, u] [A% u] - ig[A% u]d^u + g^uA^uA^^ 

- ^^A^Al - ^^u?AlAl - d^uAl - igd^u{A% a;] + d^ud^u + z| [u;, S^a;]^} + 0(0;^) 

= fA[h]-Trl d'x {a';,, d^co] + TrJd'x (/^c^o; - ^^^a;^ - ^^a;^^ 

- g'[Ay][A';^M + 5'a;A>4 - ^AyA';^ - ^^u^A^^ +Tr J d'x {d.ud.u 
+ il[u,di,u;]Ayigd^u[A>;,,u]-ig[A';,,u;]d^u + z^Ay,d^u;]) +0{u^) 
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.2 ,h 



fA[h]+2 j d^xtr {uod^A^}j + j x tr ^2g^ uj A^^uj A^^ - 2g'^ A^^A 

[a^u^ - (a> - a;^) ] + j d^x tr {d^cod^uj + i^-u^d^uAl - i^d^ucoA 



- igd^uA'l^LO + igd^.totoA''^ - igA'^^ud^u + igojA^^d^uj + i^A^uod^u: - i^A'^^d^uju^j + 0{lo^ 



fA[h] + 2TVy" d^x +TvJd^x (d^codf^u; + iguod^ujA^^ - igd^uoojAl 

2igd^LoAloj + 2igd^u:u:A^^ + 0{to^) . (A.13) 



Thus 



fA[v] = fA[h] +2Ti j d^x (c^a^4) + Tr y d^x (^d^iod^u + igud^uA'^, - igd^uJUjA'l, 
- ig (d^Lo) A^uj + ig (5^cu) c^^) + 0{ij^) 

= fA[h] +2Ti j d^x (c^a^4) + Tr y d^x {o^lo (^O^lo - ig [4, a;] ) } + 0{u;^) . 

(A.14) 



Finally 



so that 



fA[v] = fA[h] + 2Tr j d^x (c^a^4) - Tr ^ d^xvod^D^{A'')uj + 0(0;^) ^ (^.15) 



5^fA[h] > ^ > . (A.16) 



We see therefore that the set of field configurations fulfilling conditions (|A.16|) . i.e. defining 
relative minima of the functional belong to the so called Gribov region fi, which is defined 

as 

^= {A^\ d^A^ = and - d^D^{A) > 0} . (A. 17) 

Let us proceed now by showing that the transversality condition, d^A^ = 0, can be solved for 
h = h{A) as a power series in A^. We start from 

A''^ = h^A^h + -h^d,,h , (A.18) 



with 



Let us expand h in powers of 



From equation HA.18|) we have 



h = = e^^'^"^" . (A. 19) 



h = l + ig(t>- + 0(,/.=^) . (A.20) 



2 2 

4 = + ig[A^, 4>\ + 9^(t>A^<t> - y^M'A' - y</''^/. - d^ct> + i|[</., d^] + 0{<ty') . (A.21) 
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Thus, condition duA^ = 0, gives 



d^A + ig[df,Af„ </>] + ig[A^, df^cj)] + g^d^(])A^(t) + g'^(^d^A^<]) + g^<t)A^d^ 



- 3ld^A^<l?-^^A^d,^ 



g 



_ r 



9\2 



df^A^ 
(A.22) 



This equation can be solved iteratively for as a power series in An, namely 



^2^(1^^ + 1 Q2 



, ■ 9 



* g 

H — 

252 



9^ r^A 



(A.23) 



so that 



A' 



1 9 



Ay, du 



dA 

-Q2 



+ ig 



+ ^2 



1 d 

Q2 Q2 



"'292 

+ 0{A^ 



dA, ^dA 



Expression HA.24() can be written in a more useful way, given in ea. (|3.2|) . In fact 



(A.24) 



Ah 



An - ig 



Q2 



■^dA,A^ 

52 '52 



+1 



Q2 

1 

Q2' 



^_g 

2 

1 



^dA,dyj^dA 



+ 0{A^ 



-dA,d,^dA 



d d 
-^dA + ig-^dy 



-^9A, Ay 



A^ - ^dA + ig 



+ 0{A^ 
Af,,^dA 



.gdfj_ 

'202 



dA 
9^' 



dA 



Q2 

+ 0{A^ 



^g 

2 



^dA,df,j^dA 



+ ^9^2 



-^dA, Ap 



(A.25) 



which is precisely expression ()A.24|) . The transverse field given in eq. (|8.2|) enjoys the prop- 
erty of being gauge invariant order by order in the coupling constant g. Let us work out the 
transformation properties of (py under a gauge transformation 



5A^ = -d^LJ + ig[Af„uj] 



(A.26) 



We have, up to the order 0{g^ 



-dyuo + ig 

-dyUJ + l- 



-^dA,dyUJ 
^dA,dyUJ 



u),dy-^dA 



■g 



dA 

-Q2 



,dyUJ 



+ Oig' 



+ ^2 



dy-^dA,Lo 



Therefore 



-dy 



UJ — l- 



dA 
92"' 



UJ 



+ 0{g^). 

+ 0{g''). 



from which the gauge invariance of A^ is established. 



(A.27) 
(A.28) 
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Finally, let us work out the expression of A'^- as a power series in An. 



Tr / 



Tr / df^x 
Tr / d'^x 



A„ 



Q2 

-^(^A, A^ 



+ 



^9 



1 

92 



Q2 



Tr d'^x -lA 



Ay - ig 



1 

-^dA, Ay 



^_9 
2 



dA,d^^dA 

1 1 

-^dA,dy-^dA 



'Ml' 



A,^ 



( OA 
2ig lAy- dy-^ 



OA 



■ ( . ^ dA 
+ig ( Ay - Oy-^ 



Tr d^x -lA 



dA dA 
df,dy 



+ 0{A' 



+igAy 



dA „ dA 
92 



Q2 

dydA 



Ay - 2igAy 



dA 



+ 2ig- 



dydA 



Q2 



dA ■ 

-^^Ay 



92 



dA dA 



+ 0{A' 



d^x 



d^,dy 

Q2 



92 52 



^2 



leading to the result quoted in eq. (|3.4|) . 



+ 0(A) . 
(A.29) 



We conclude this Appendix by noting that, due to gauge invar iance, ^^jjj can be rewritten 
in a manifestly invariant way in terms of F^^y and the covariant derivative Df^ [nU] , see eq. ()3.17p . 



B Properties of the Stueckelberg term. 

In this Appendix we derive some useful properties of the non-Abelian Stueckelberg term 0$ 
p] . defined by the equations (|SIII)-(|SIH1)- The expression (|3.7|) is left invariant by the gauge 
transformations given in eq. (|3.9|) . In fact 

(^A^-'-U-^d^U^ ^ V-^(^A^--U-^d^U^V . (B.l) 

Thus 

Ti{A^--U-^d^U^ ^ Ti(^A^--U~^d^U^ . (B.2) 

Let us look now at the equations of motion of the Stueckelberg field 0", as expressed in eq. (|3.11)) . 
from which 

d^A^ - -d^ {u-'d^u) - [A^, u-^d^u] = o . (b.3) 
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Expanding the term U ^dnll in power series of 



2 2 

igT^d^cp" - y {d^cp") cP^T^T^ - ^cP^d^^cp^T^'T^ + g'^cP'^d^cp^T'^T'' + 



= igT'^d.r - y {d>.r) T'] + 0{cP^) , (B.4) 

U~'d^U = igT^d^r - ^iT'^r'^ {d^r) + 0{<p') . (B.5) 
After substitution of expression (|B.5|) in ea. (|B.3|) . we have 

= a^A^lT" + a^ (r'^d^r - ^T'^r'"'d^r<p') - igAl,d^^'[T\ T^] + higher order terms 
= T" i^dA" + d'^(p'' - |/"''" (a^/) (/.^ + gf^^'^Ald^^'''^ + higher order terms , (B.6) 



yielding 



from which 



d'^cP" = -dA'' - gf^^A^^df^cp^ + 1/'^''^ (aV^) (/^^ + higher order terms , (B.7) 



^dA- - ^ (r'^'^A%cP'^ - ^-r''^ (a^/) <p^ ) + higher order terms . (B.8) 



and 

<P'' = -^dA- - ^ i^r^A'^d^cP' - I ^-'^ '^'-'^ 
Finally, substituting recursively for (p, we obtain the expression H3.12() . 
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